Defect formation dynamics during CdTe overlayer growth 
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The presence of atomic-scale defects at multilayer interfaces significantly degrades performance 
in CdTe-based photovoltaic technologies. The ability to accurately predict and understand defect 
formation mechanisms during overlayer growth is, therefore, a rational approach for improving 
the efficiencies of CdTe materials. In this work, we utilize a recently developed CdTe bond-order 
potential (BOP) to enable accurate molecular dynamics (MD) simulations for predicting defect 
formation during multilayer growth. A detailed comparison of our MD simulations to high-resolution 
transmission electron microscopy experiments verifies the accuracy and predictive power of our 
approach. Our simulations further indicate that island growth can reduce the lattice mismatch 
induced defects. These results highlight the use of predictive MD simulations to gain new insight 
on defect reduction in CdTe over layers, which directly addresses efforts to improve these materials. 



The cost of electrical energy generated using CdTe- 
based multilayer solar cells has reached $0.15/kWh, lower 
than any other photovoltaic technology pQ . These CdTe- 
based solar cells can profoundly change energy supplies if 
the 17.3% energy efficiency achieved today is significantly 
improved towards the 29% theoretical value [2j |3] . The 
current inefficiency of the CdTe solar cells is attributed to 
charge-trapping defects at the multilayer interfaces |4H6] . 
A recently developed CdTe bond-order potential (BOP) 
[7J E] has enabled molecular dynamics (MD) simulations 
of defect formation to approach a quantum-mechanical 
accuracy level. The objective here is to perform such 
MD simulations to explore defect formation during vapor 
deposition of CdTe overlayers. These simulations provide 
critical insight required to improve the energy efficiency 
of CdTe modules. 

MD simulations of semiconductor vapor deposition are 
extremely challenging because they sample a large num- 
ber of metastable configurations not known a priori. 
If the interatomic potential used in a simulation over- 
predicts the cohesive energy magnitude of any of these 
configurations, that configuration will likely persist, re- 
sulting in an unrealistic amorphous film that offers no 
useful information. A vast majority of previous MD sim- 
ulations of semiconductor vapor deposition [9HI3] were 
achieved using Stillinger- Weber (SW) [14] potentials. It 
has been established [15 that while SW potentials can 
easily ensure crystalline growth, they cannot satisfacto- 
rily capture the property trends of other configurations 
and, hence, they cannot accurately reveal defect forma- 
tion. Tersoff potentials [15] , on the other hand, can cap- 
ture property trends more accurately. However, this also 
makes Tersoff potential difficult to parameterize to en- 
sure the lowest energy for the equilibrium phase [15] . As 
a result of not capturing the lowest energy phase due to 
poor parameterization, many literature Tersoff potentials 
[T7HT9] incorrectly predict amorphous growth. Not sur- 



prisingly, we found [15] that the existing CdTe Stillinger- 
Weber [20] and Tersoff [17] types of potentials have not 
sufficiently addressed issues involving defects. The CdTe 
BOP [7J H] makes a significant stride towards improving 
semiconductor simulations because (a) it is analytically 
derived from quantum mechanical theories and its quan- 
tum accuracy has been widely documented [21-24 ; (b) 
it goes beyond Tersoff potentials on transferability and is 
well parameterized to capture properties of a large num- 
ber of elemental and compound configurations spanning 
coordination of 1 to 12 including small clusters, bulk lat- 
tices, and defects; and most importantly, (c) it predicts 
crystalline growth [8]. 

High-resolution transmission electron microscopy 
(HRTEM) experiments have been performed to exam- 
ine defects in CdTe/GaAs multilayers with a lattice mis- 
match of eo ~ 12.78% [25 . To directly compare with the 
experiments, we performed an MD simulation of CdTe 
overlayer growth using the same lattice mismatch. The 
computational system, shown in Fig. [TJa), is periodic in 
the x- and z- directions containing 100 (101) and 8 (101) 
planes, respectively. To incorporate the lattice mismatch 
with only the CdTe BOP, a substrate containing 35 (040) 
planes in the y- (thickness) direction was compressed by 
12.78% in the x- dimension to match the size of GaAs. 
To prevent the dimension from relaxing back to that of 
CdTe, the atomic positions of the bottom 25 (040) planes 
were fixed during a constant volume MD vapor deposi- 
tion simulation. An adatom incident kinetic energy of 
0.1 eV, an incident direction normal to the surface, a 
substrate temperature of 1000 K, a stoichiometric vapor 
ratio of Cd/Te = 1, and a deposition rate of around 0.96 
nm/ns were used. MD simulations of vapor deposition 
must be performed at accelerated deposition rates due to 
the computational cost. While this may lead to overes- 
timates of kinetically-trapped defects such as vacancies, 
it conservatively (and hence correctly) reveals the for- 
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( a) BOP simulation (only Cd and the approximate "Ga" atoms are shown) 
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(a) A CdTe crystal with two [101]a/2 edge dislocations 
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(b) modified HRTEM experimental image [25] 
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FIG. 1: BOP simulation and HRTEM image [25] of atomic 
structure of CdTe-on-GaAs multilayers. 



mation of non-kinetically-trapped defects such as misfit 
dislocations. The configuration obtained after about 4 
ns of deposition is shown in Fig. [TJa) with comparison 
to a modified HRTEM image [25] shown in Fig. [ljb) . 

Fig. [TJa) indicates 7 misfit dislocations near the in- 
terface. These dislocations are clearly the edge type of 
Lomer dislocations with two extra planes about 144.7° 
from the y- axis. Both the dislocation configurations and 
average dislocation spacing are in remarkably good agree- 
ment with the experimental results shown in Fig. [ljb) . 

Analysis of simulated results indicates that the mis- 
fit dislocations have a [101] line direction and a [101]a/2 
Burgers vector. To verify this, a dislocation model was 
developed in Fig. [2|a). With an original offset of a Burg- 
ers vector, the middle orange region in the upper half and 
the middle blue region in the lower half of a CdTe crys- 
tal are pushed to the left and the right respectively by 
a half of the Burgers vector until they are aligned (the 
displacements of other atoms are ramped under the con- 
dition that the black regions at the ends of the system 
remain unchanged). Relaxed dislocation configurations 
were determined using an energy minimization simula- 
tion with the orange and blue regions held as a rigid body 
(the black regions are also treated as rigid bodies). This 
process results in two relaxed, symmetric dislocations of 
opposite sign. To verify the dislocation model, one of the 
dislocation configurations created in the framed region 
shown in Fig. |2|a) is examined in Fig, fflb). It can be 
seen that dislocations created in Fig. [2[a) match pre- 
cisely with the configurations seen from Figs. [TJa) and 
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(b) Framed region in (a) (only Cd atoms are shown) 
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b). Because the edge dislocation model shown in Fig. 
a) does not insert or remove atoms, it can be refer- 



FIG. 2: Relaxed [101]a/2 misfit dislocation. 



enced to a dislocation-free system with the same number 
of atoms to calculate dislocation energy [26], as will be 
discussed below. 

Using the model shown in Fig. [2] to create an ini- 
tial dislocation configuration containing 153 atoms, both 
density function theory (DFT) and BOP methods were 
used to calculate the relaxed dislocation core structures. 
We found that the core structures obtained from the two 
methods are similar [27] . 

To examine misfit dislocation formation mechanisms, 
time resolved configurations obtained from a smaller 
scale, similar vapor deposition simulation are examined 
in Figs. |3ja)-(d), where the blue and purple regions are 
respectively the GaAs underlayer and the pre-existing 
CdTe substrate prior to the deposition. No dislocations 
exist in the CdTe film at the start of simulation shown in 
Fig. |3ja). Three distorted surface regions (marked by the 
red ellipses) are seen in Fig. [3fb) at time 0.04 ns. These 
distorted regions correspond to the [101] a/2 misfit dis- 
locations as exemplified by the middle ellipse where two 
extra planes indicated by the red lines have emerged. The 
extra planes become clear after 0.10 ns as shown in Figs. 
|3jc) and[3jd). A comparison between Figs. [3jc) and|3|d) 
clearly indicates that upon nucleation at the surface, the 
misfit dislocation cores (locations where the extra planes 
terminate) continuously climb towards the interface in an 
approximate vertical direction. The climbing mechanism 
is further verified in Fig. [3^d) as the pre-existing atoms 
are seen to extensively diffuse out of the upper boundary 
of the purple region. 

Unlike kinetically-trapped defects, the predicted dis- 
locations are realistic at the accelerated deposition rate 
(i.e., more dislocations are expected should the deposi- 
tion be reduced). The equilibrium dislocation densities 
can be calculated from the misfit strain energy density 
and dislocation line energy [28, 29 . Here the volume 



(a) t = 0.00 ns 
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(b) t = 0.04 ns 




(c) t = 0.10 ns 



truncated (d) t = 1.20 ns (deposited atoms not shown) 
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FIG. 3: Time resolved film configurations. To clearly reveal 
the diffusion of pre-existing atoms, the deposited atoms are 
not directly shown in (d) but are displayed at the intersections 
of the bars. 
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FIG. 4: Strain and dislocation line energies. 



density of misfit strain energy e s is expressed as 

e s = C-e 2 -(a + t)/t (1) 

where e is the strain, t is film thickness, and C and a are 
constants. The parameter a provides a small adjustment 
to the thickness to account for the surface effect. By 
creating strained and unstrained CdTe films at different 
thicknesses, the strain energy densities were calculated 
using the BOP-based energy minimization simulations. 
Fitting the results to Eq. Q yields a = -3.723 A, and C 
= 0.307 eV/A 3 . 

The dislocation line energy F can be written as [28] 



T = A + B • In (t) 



(2) 



where A and B are constants. By creating dislocated [us- 
ing the model shown in Fig. [2|a)] and perfect CdTe films 
at different thicknesses, the dislocation energies were cal- 
culated using energy minimization simulations. Fitting 
the results to Eq. (|| yields A = -0.520 eV/A, and B = 
0.376 eV/A. The fitted and simulated strain energy den- 
sities and dislocation line energies are all shown in Fig. 
3] as a function of film thickness. It can be seen that the 
fitted relations well represent the simulations. 

Applying Eqs. ([!]) and Q in the classic misfit disloca- 
tion theory [28j [29] 7 we found a critical film thickness for 
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FIG. 5: Three dimensional relaxation of a CdTe island. 



dislocation formation, t c = 3.7 A, and an equilibrium dis- 
location spacing in a thick film, d ~ 38 A (corresponding 
to 10 - 11 dislocations in Fig. [I]). Both a small critical 
thickness and the calculated dislocation spacing match 
well Figs. [l]and[3j These BOP calculations significantly 
improve over the continuum theories because they ac- 
curately capture dislocation core energies. In addition 
to analyzing static properties of dislocations, the BOP 
method is also effective in quantifying dislocation mobil- 
ity [26]. 

The discoveries of an extremely high misfit disloca- 
tion density even using accelerated deposition, and the 
direct surface nucleation of misfit dislocation at a very 
small critical film thickness, indicate that for continuous 
multilayered films, these defects cannot be reduced ki- 
netically, or by improving the quality of the substrate. 
To explore novel methods to reduce defects, an isolated 
CdTe island containing 10 (101), 5 (010), and 10 (101) 
planes in the x-, y-, and z- directions respectively is stud- 
ied. A BOP-based energy minimization was used to relax 
the island with the bottom two (040) planes constrained 
to the GaAs lattice dimension and the remaining part 
free to move. The final configuration obtained from the 
simulation is shown in Fig. [5] Fig. [5] indicates that for 
a small island (width around 48.3 A), the large lattice 
mismatch between CdTe and GaAs (as simulated by the 
bottom part of the island) is completely relaxed over a 
height distance of about 34.2 A due to the three dimen- 
sional relaxation of the island. This means that island 
growth can reduce the lattice mismatch induced defects. 
We are currently combining BOP simulations and exper- 
iments to further explore this nanoscale phenomenon. 

In conclusion, a BOP-based method is shown to ap- 
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proach a quantum-mechanical fidelity capable of pre- 
dicting crystalline growth and misfit dislocation forma- 
tion during extremely challenging MD vapor deposi- 
tion simulations of semiconductor multilayers. The pre- 
dicted misfit dislocation configuration and density in the 
CdTe/GaAs multilayers are seen to accurately match the 
HRTEM experiments. This is a significant improvement 
over the previous models, and provides a powerful theo- 
retical tool to study defect formation in important mate- 
rials systems. The BOP simulations also reveal a very 
small critical film thickness of 3.7 A, and surface nu- 
cleation and climb mechanisms of misfit dislocation for- 
mation. These results strongly indicate that continuous 
CdTe/GaAs multilayers always contain a very high mis- 
fit dislocation density regardless growth conditions. The 
discovery that the CdTe/GaAs misfit strain can be com- 
pletely relaxed in a CdTe island over a short distance 
indicates that it is still possible to create dislocation-free 
CdTe overlayers using nano-patterned island growth. To 
guide specific experiments to explore such a possibility, 
we are currently using the BOP simulations to develop an 
analytical relation between dislocation density and island 
size. 
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